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Abstract. We address realistic schemes for the generation of non-Gaussian 
states of light based on conditional intensity measurements performed on 
correlated bipartite states. We consider both quantum and classically correlated 
states and different kind of detection, comparing the resulting non Gaussianity 
parameters upon varying the input energy and the detection efficiency. We find 
that quantum correlations generally lead to higher non Gaussianity, at least in 
the low energy regime. An experimental implementation feasible with current 
technology is also suggested. 



1. Introduction 

Continuous variable (CV) quantum information has been developed with Gaussian 
states and operations, which allow the realization fundamental protocols as 
teleportation, cloning and dense coding, and still play a relevant role in CV quantum 
information processing. Recently, also the non Gaussian (nG) sector of the Hilbert 
space is receiving attention for the potential application in entanglement distillation 
and entanglement swapping for long distance quantum communication. In turn, using 
nG states and operations, teleportation [TJ[2l[3] and cloning [4] of quantum states may 
be improved. In order to quantify the amount of nG of a state measures have been 
recently proposed, based on distances between the quantum state under investigation 
and a reference Gaussian state [7] . 

Gaussian states are generated using linear and bilinear interaction in optical 
materials, whereas the generation of nG states requires higher nonlincarities as those 
exhibited by a Kerr medium. An alternative approach is to exploit the effective 
nonlinearity induced by conditional measurements. In fact, if a measurement is 
performed on a portion of a composite system, the other component is conditionally 
prepared according to the outcome of the measurement and the resulting dynamics 
may be highly nonlinear. The the rate of success in getting a certain state is 
equal to the probability of obtaining a certain outcome and may be higher than 
nonlinear efficiency, thus making conditional schemes possibly convenient even when 
a corresponding Hamiltonian process exists. 
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In this paper we focus on preparation schemes feasible with current technology 
and address generation of non-Gaussian states of light by conditional intensity 
measurements on the twin-beam state (TWB) generated by parametric down- 
conversion and on classically correlated states generated by mixing a thermal beam 
with the vacuum in a beam splitter (BS). In particular we consider different kinds 
of conditional measurements: ideal photodetcction, inconclusive photodctcction and 
photodctcction performed with a finite quantum efficiency rj and no dark counts. We 
compare the resulting nG varying all the involved parameters and in particular the 
number of photons of the initial state. 

Let us consider a CV system made of d bosonic modes described by the mode 
operators afc, k = 1 . . . d, with commutation relations [afc,aj] = Skj- A quantum 
state q of d bosonic modes is fully described by its characteristic function x[g](A) = 
Tr[p -D(A)] where D{X) = <S)^ =1 Dk{\k) is the d-mode displacement operator, with 
A = (Ai, . . . , Xd) T , Afc e (D, and where Dk(Xk) = exp{Afcdj. — \* k ak\ is the single- 
mode displacement operator. The canonical operators are given by q k = (fl/t + a\)/\/2 
and pk = (at — al.)/\/2i with commutation relations given by [qj,pk] = iSjk- Upon 
introducing the vector R = (q±,px,.. . ,qd,Pd) T , we have the vector of mean values 
X = X[g], Xj = (Rj), and and the covariancc matrix (CM) er = er[g] of a quantum 
state, o k j = \{RkRj + R 3 Rk) - (Rj}{Rk) where (O) = Tr[g O]. A quantum 
state qg is said to be Gaussian if its characteristic function is Gaussian, that is 



X[£>g](A) = cxp | -iA T erA + X T Al, where A = (RcAi,ImAi, . . . , ReA d , lm\ d ) T . A 



Gaussian state is fully determined by CM and X. 

For a generic CV quantum state g, a measure of nG based on the quantum relative 
entropy has been introduced in [6] as 5[g] = S(t) — S(g) where S(g) = — Ti[g\ogg] 
is the von Neumann entropy of a quantum state g, and r is the Gaussian state 
with the same CM and X as g. In this paper we will deal with single-mode 
non-Gaussian states that can be written as diagonal mixtures of Fock states g = 
12^=0 Pn\ n )( n \i f° r which the reference Gaussian state is the single-mode thermal state 
v N = 1/(1 + N ph ) J2n=o ( N ph/(1 + N ph)) n \n) (n\ where N ph = Tr[g a^a] is the mean 
photon number. The von Neumann entropy for these states can be easily obtained 
and therefore the nG can be written as 



2. Non-Gaussian states by conditional measurements 

We want to study the quantum states generated by conditional measurements on 
quantum and classically correlated states. In particular as initial states we consider 
(i) an entangled TWB, obtained from (spontaneous) parametric down-conversion 
(SPDC) in second-order nonlinear crystals, (ii) a thermal state mixed with the 
vacuum state in a BS of transmissivity T. The TWB state can be written as 
|A)) = s/1— |A| 2 J2 n X n \n)i (g> \n)2 where \n)j denotes the Fock number state in 
the Hilbert space of the jth mode. The parameter |A| < 1 may be taken as real 
without loss of generality and the mean number of photons for each mode is given by 
N = Tr[|A))((A| ni g> 1] = Tr[|A))((A| l<gm 2 ] = A 2 /(l - A 2 ) The output state from a 
BS of transmissivity T mixing a thermal state z^jv and a vacuum, is described by a 




(1) 
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density operator 

R = £ T (^)/ 2(1 _ T)(^)/ 2 , p+s , t+ ^^ + s ) (*^)l*><*l ® |P><g| (2) 

where i//^ = <5, l;fc (l + 2A0~ 2 [2A7(1 + 27V)] fc and 2iV is the mean photon number in R. 
In the next subsection we will derive the states generated from \X}} and R by means 
of different conditional measurements, that is 

Tr 2 [|A))((A| 1«M(')] Tr 2 [-R 1 <g> AfW] 



SA,i = m[ |,,,,,,| 1 , ~ : tuy\ 0R 



Tr[|A))((A| 1 ® AfW] Tr[R 1 <g> ¥<•>] 

where i denotes the kind of measurement we are performing and AfW is the operator 
of the corresponding probability operator valued measure (POVM). In particular we 
will consider (1) ideal photodetection, (2) (ideal) inconclusive photodetection and (3) 
photodetection with a non-unit quantum efficiency. We will then evaluate their non-G 
parameter and describe their photon statistics. In particular we will give the mean 
photon number N p h, the variance a^ h and the corresponding Fano Factor defined as 

2.1. Ideal photodetection 

The POVM of an ideal photon-resolving detector is given by projectors on the Fock 
number basis: P m = \m){m\. By taking the TWB |A)) as the initial state, the output 
state after measuring m photons is simply the corresponding Fock state £a,i = \m)(m\. 
The photon statistics of £a.i is trivial, since N p h = m and both a^ h and F p h are 
null. The nG of the output state can be easily evaluated by using Eq. ([1]) and it 
monotonically depends on the number of detected photons m. 

In the case of the classically-correlated state R with N mean photons in each 
mode, performing an ideal photodetection on one mode leaves the other mode in the 
following conditional state: 



QR.l 



l + 2N(l-T) \ m+1 ^fm + s\ ( 2NT ] 

1 + 2N M m (iT^I |S)(S| 



In this case the mean number of photons, the variance and the Fano factor of the 
ouput state are non-trivial but can be evaluated obtaining the following results 

2(1 + m)NT 2 _ 2{1 + m) NT (1 + 2N) 

ph ~ 1 + 2JV(1-T) aph ~ (1 + 2N{1 -T)) 2 



1 + 2N 



ph 



1 + 2JV(1-T) 

The corresponding nG #[f?.R,i] can be evaluated numerically and is plotted in Fig. [1] 
along with the nG of the previous case. As one can observe from the plot, the state 
generated by the classically-correlated state is less non-Gaussian than that generated 
by TWB. Its nG increases with the number of detected photons m, with the initial 
mean number of photons N and with the transmissivity of the BS T. Though 
increasing with N, we have numerical evidences that the asymptotic value for N — > oo, 
is below the nG of the Fock states \m)(m\ obtained from the TWB. 



Non- Gaussian states by conditional measurements 



4 



5[P] 



1.5_ 

1.0 

0.5 



1 .0 \ 
0.8 



6 8 10 



Figure 1. Left: Non-Gaussianity of quantum states generated by performing an 
ideal photodetection on one arm of a TWB A}) (dashed lines) or of the classically 
correlated state R (solid lines) obtained with a BS with transmissivity T = 0.9, 
as a function of the initial average number of photons N and for different values 
of detected photons m. From lighter to darker gray: m = {1,2}. Right: Non- 
Gaussianity of quantum states generated by inconclusive photo-detection on a 
TWB |A)) (black dashed line) or on the quantum state R (solid lines) for different 
values of the BS transmissivity T (from lighter to dark gray: T = {0.5, 0.99}) as 
a function of the initial number of photons N . 



2.2. Inconclusive photodetection 

The ideal inconclusive photodetection is described by the POVM operator H ff = 
|0)(0|, if no photons are detected, and H on = 1 — |0)(0| if one or more photons are 
detected. Let us first consider the state obtained when one ore more photons are 
detected on one mode of the TWB |A)). The state that is generated can be written 
in the Fock basis as 

s=l 

Note that, for A — > (TV — ► 0) we obtain the Fock state while for A — * 1 

(TV — ► oo), the state is not physical. The photon statistics of qr,2 can be summarized 
in terms of the initial number of photons TV using the parameters N p h = 1 + N, 
a ph = N(l + TV), and F p h = TV. The corresponding nG <S[gA,2] can be evaluated 
numerically and is plotted in Fig. [TJ It is apparent from the plot that 5[f?A,2] is a 
monotonically decreasing function of the initial number of photons TV, starting from 
the nG of the Fock state |1)(1| and approaching zero for TV — ► oo. 

Let us consider now the conditional state generated from the classically-correlated 
state R. It is diagonal in the Fock basis and can be written as 



1 + 27V(1-T) 

QR,2 - 



( 2NT ^ 


S 




( 1 + 27V > 


*- 1 1 


V 1 + 2NJ 


1 + 2NT 


Vi + 2A/ry 


l + 2iV_ 



27V (1 - T) 

The average photon number, its variance and the Fano factor are 



\s)(s\ 



47VT(1 + 7V(1-T)) 2 _ 47VT(1 + 7V(3-T + 4T7V 2 (1 - T) 2 + 27V(1 - T 2 ))) 

ph ~ l + 2/V(l-T) ° vh ~ (1 + 27V(1-T)) 2 

Fph ~ 2(1 + N T) 2(1 + 7V) 1 + 27V — 

p y ' 1 + 7V(1 - T) 1 + 27V(1 - T) 

The nG 5[qr.2\ can be numerically evaluated and is plotted in Fig. [T] In this case it 
has not a monotonic behavior as a function of TV while it is still an increasing function 
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of the transmissivity T. In the low energy regime the state generated starting with 
a TWB is definitely more nG than that generated by means of R. By increasing the 
initial mean number of photons, since approaches zero, we observe a region 

where the classical correlated state gives birth to a more non-Gaussian state than the 
quantum one. 



2.3. Inefficient photodetection 

The POVM of a photon-number resolving detector with quantum efficiency rj and no 
dark counts is given by the Bernoullian convolution of the ideal number projectors 
Pi = \l)(l\, and thus, for m detected photons, it is described by the operator 
n m = V m J2hL m (l ~ 7 l) l ~ m (m)P l - The state generated by detecting m photons on 
one arm of a TWB can be thus written as 



£A,3 



l-A 2 (l-7y) 
(A 2 (l - 77))" 1 



E 



(\ 2 (l-r,)) l \l)(l\. 



Here we summarize the photon statistics of this state in terms of the previously 
introduced quantities, by substituting A with N 



(l+m)(l + AQ 

Mph = : ~ t 



1 + Nrj 



Fph — 



'ph 



(l + m)(l -r])N(l + N) 



(1 + vNy 



N(l + N)(l + m)(l-r)) 



(m(l + N) + N{1 - rf))(l + Nrf) 

The nG can be evaluated numerically and it is plotted in Fig. [5] for different values 
of the parameter. We observe that, as expected, nG is a monotonically increasing 
function of the number of detected photons m. Moreover we observe that it decreases 
with the initial mean photon number N, while it increases at increasing detector 
efficiency. By choosing 77 = 1 we trivially obtain the results of the ideal photodetection. 
For N — > 00, the state generated by measuring a TWB turns out to be the 
corresponding (normalized) POVM operator, that is £»a.3 = rill m . The asymptotic 
values approached for N — *• 00 correspond to the nG of this particular state, and 
increases at increasing r\. Finally we consider the classically-correlated state R. Also 




Figure 2. nG of conditional states generated by performing an inefficient 
photodetection on one mode of a TWB |A)) (left figure) or of the classical 
correlated state R (right figure) obtained by a BS with transmissivity T = 0.5, 
as a function of the initial average number of photons iV, for m = 2 photons 
detected and for different values of the efficiency 77. From darker to lighcr gray: 
7] = {0.8, 0.9, 0.99}. As one can observe comparing the two plots, the nG of states 
generated from TWB are higher than the one obtained from thermal states. 
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in this case we can obtain the state generated by measuring m photons in one arm, 
with efficiency r\ 

( l + WnO-T) \ m ^^fm + s\ ( 2NT V 

QR.,3 y l + 2N ( T + r] (i- T ))J m J \l + 2N(T + r,(l-T))J l ' V[Sl 

The mean number of photons, its variance and Fano factor result 

2(1 + m)NT 2 _ 2(1 + m)NT{l + 2N(T + - T))) 

ph ~ 1 + 2N V (1 - T) aph ~ (1 + 2Nn(l — T)) 2 

2NT 



i 



1 + 2Nr}{\ - T) ' 

Again the nG <5[gfl,3] has been numerically evaluated and is plotted in Fig. O As in 
the ideal case, the nG obtained from a classically-correlated state take lower values 
than those from TWB. Moreover we again observe that a higher nG is obtained at 
higher transmissivity (as in the previous cases) and, unlike the TWB case, at lower 
detection efficiency and higher values of the initial number of photons N . For N — > oo 
the asymptotic state can be written as 



QR,3 



r?(l-T) \ m+1 ^fm + s\f T 



T + V (1-T)J ^\ m J\T + r){l-T) 



Again we have numerical evidences that for N — > oo the nG obtained for the classical 
correlated state are below those obtained by measuring the TWB state. 



3. Experimental proposal 

The experimental implementation of the optical states described in the previous 
sections, namely TWB states and thermal states, would allow us to verify the 
correctness of our model. Obviously, as the real detectors are endowed with a non-ideal 
quantum efficiency, here we consider only the case of imperfect detection presented in 
Sec. [2] In order to match the requirements of the theoretical model we plan to use 
two hybrid photodetectors (HPD, Hamamatsu) as the detectors, which are endowed 
with a partial photon- number resolving capability and no dark counts [10| . 

First of all, we discuss the case of the generation of TWB states with a sizeable 
mean photon number. Such states can be obtained by pumping a second-order 
nonlinear crystal (cither in typc-I or typc-II phase-matching) by means of a pulsed laser 
in a travelling- way configuration [8] . As sketched in the upper panel of Fig. [3l we can 
exploit either the third harmonics (A = 349 nm, 4.4-ps-pulse duration) or the fourth 
harmonics (@ A = 262 nm, 4-ps-pulse duration) of a frequency-tripled Nd:YLF mode- 
locked laser amplified at 500 Hz to produce bright cones of spontaneous parametric 
down-conversion in a type-I /3-BaB204 crystal. The TWB state will be measured 
by selecting with suitable pin-holes a pair of twin coherence areas on the signal and 
idler cones [|5]. In order to make the quantum efficiencies of the two arms as equal 
as possible, we will operate at frequency degeneracy. Note that with both choices 
for the pump, the measured TWB wavelength will be in the visible spectrum range, 
where most photodetectors have the maximum quantum efficiency. Multi-mode fibers 
having a good transmissivity in the visible spectrum range will be employed to collect 
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Figure 3. Proposed experimental setups. Upper panel: TWB state. Lower 
panel: classically correlated thermal state. 



all the light passing the pin-holes and to send it to the detectors. The study of the 
dependence of the nG factor on the overall quantum efficiency will be experimentally 
reproduced by inserting a variable neutral density filter on the arm of the detector 
performing conditional measurement. Unfortunately, depending on crystal length and 
phase matching conditions, the TWB state produced by ps-pulsed lasers is intrinsically 
temporal multi-mode. For this reason, the single-mode theory previously developed 
should be extended to correctly describe the experimental situation. 

The production of a single-mode pseudo-thermal state can be easily obtained by 
inserting a rotating ground glass plate on the pathway of a coherent field, followed by 
a pin-hole selecting a single coherence area in the far-field speckle pattern [11] (see the 
lower panel of Fig. [3]). To match the maximum quantum efficiency of the detectors, 
we can use the second harmonics (A = 523 nm, 5.4-ps-pulse duration) of the Nd:YLF 
laser described above. The thermal light will then be split in two parts by means of a 
BS in order to perform conditional measurements. To make the setup more versatile 
and in particular to have the possibility to study the dependence of the nG factor 
on the transmissivity of the BS, it is convenient to substitute the ordinary BS with a 
system composed by a half-wave plate (A/2 in Fig. [3]) and a polarizing cube BS. Again 
the measurements can be performed with a pair of hybrid photodetectors. 

4. Conclusions 

We have suggested an experimentally feasible scheme for the generation of non- 
Gaussian states of light by conditional intensity measurements on correlated bipartite 
states. We have considered both quantum and classically correlated states and 
evaluated the amount of resulting nG for different kind of detection and input energy. 
Although the emergence of nG is not related to the quantum nature of the involved 
beams we found that quantum correlations generally lead to higher nG, especially in 
the low energy regime. 
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